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Abstract: This paper focuses on the design of Generic Model 
Control (GMC) algorithm for the control of the acetate 
concentration in a biomethanation process with production of 
methane gas. A nonlinear dynamic model for the 
biomethanation process is used in the derivation of the control 
law of Generic Model Control scheme. The controller 
parameters are estimated using the Trial and Error method. 
The effectiveness and performance of the proposed control 
strategy is illustrated by a MATLAB program under various 
operation conditions and disturbances. The simulation results 
obtained confirmed the good quality of the control. 
Keyword: Generic model control; process control; 
biomethanation process 

I. Introduction 

The control of biotechnological processes has been an 
important problem attracting wide attention. The main 
engineering motivation in applying advanced control methods to 
such processes is to improve operational stability and production 
efficiency. But, the use of modern control for these bioprocesses 
is still low. The nonlinearity of the bioprocesses and the 
uncertainty of kinetics impose the advanced control strategy as a 
suitable approach. So, the difficulties encountered in the 
measurement of the state variables of the bioprocesses impose 
the use of the so-called "software sensors". Note that these 
software sensors are used not only for the estimation of the 
concentrations of some components but also for the estimation of 
the kinetic parameters or even kinetic reactions. The dynamics of 
this biotechnological process are described by a set of nonlinear 
differential equations obtained from the reaction scheme and the 
unknown reaction rates are estimated. For the estimation of 
unknown parameters of the process, the distribution approach is 
used. The parameter identification of deterministic nonlinear 
continuous -time systems (NCTS), modeled by polynomial type 
differential equation has been considered by numerous authors 
[1], [2]. The generic model control scheme used in this paper is 
based on a control law that uses the nonlinear dynamic model of 
the process to calculate the effect of sequences of control steps 
on the controlled variables. 

The paper is organized as follows. Section II is devoted to 
description and modeling of a biomethanation process. The 
generic model control strategy is presented in Section III. The 
generic model control law design is derived in Section IV. 
Simulations results presented in Section V illustrate the 
performance of the proposed control algorithms and, finally, 
Section VI conclusion of the paper. 



II. Modeling of a Biomethanation Process Description 

Anaerobic digestion is a multi-stage biological wastewater 
treatment process whereby bacteria, in the absence of oxygen, 
decompose organic matter to carbon dioxide, methane and water. 

In the first phase, the glucose from the wastewater is 
decomposed in fat volatile acids (acetates, propionic acid), 
hydrogen and the second phase; the ionized hydrogen 
decomposes the propionic acid CH3CH2COOH in acetates, H2 
and carbon dioxide C02. In the first methanogenic phase, the 
acetate is transformed into methane and C02, and finally, in the 
second methanogenic phase, the methane gas CH4 is obtained 
from H2 and C02, [5]. The following simplified reaction scheme 
is considered [1], 



01 

51 > x x + s 2 

01 

5 2 > X 2 + P 



In this paper, we consider a biomethanation process 
wastewater biodegradation with production of methane gas that 
takes place inside a Continuous Stirred Tank Bioreactor whose 
reduced model is presented in [1], [2]. 

The dynamic equations of the biomethanation process 
wastewater biodegradation are derived using material balance 
relations regarding, Si represents the glucose substrate, S 2 the 
acetate substrate, X 1 is the acidogenic bacteria, X 2 the 
acetoclastic methanogenic bacteria and P represents the product, 
i.e. the methane gas. The reaction rates are denoted by 0i , 0 2 
and Q is the methane gas outflow rate. D represents the dilution 
rate and S; represents the concentration of the influent substrate - 
glucose. This process is evidently highly non-linear due to the 
inter-relationships of the states [1], [2]. According to the general 
material balance equation, the mass balance equations are 
represented by equations (1-5). 



dX 1 

^=01-^1 

-±= -K 1 Q 1 -DS 1 +DS in 
dX 2 

~^ = 0 2 -DX 2 
dS 2 

-^ = K 2 0 1 -K^ 2 -DS 2 
dP 

— = K 4 0 2 -DP-Q 



(1) 
(2) 
(3) 
(4) 
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III. Strategy of Generic Model Control (GMC) 

GMC is an advanced model-based control strategy which 
uses linear/nonlinear models of a system to compute the control 
action. Since the GMC can directly use nonlinear models of a 
process to determine a control action, the nonlinear models do 
not need linearization. Generally the GMC controller design is 
based on the following nonlinear state space [3-5] 

dx 

= x =f(x) + g(x)u (6) 



dt 



y = h(x) 



(7) 



Where /(x) and^f(x) are vector fields, i.e. they are vector 
valued functions of a vector and h(x) is a scalar field, i.e. a 
scalar valued function of a vector x . Generic Model Control 
(GMC) uses a model of the process in formulating the control 
law. However, rather than adopting a classical approach of 
comparing the trajectory of the process output against a desired 
trajectory, GMC defines the performance objective in terms of 
the time derivatives of the process output, i.e. minimizing the 
difference between the desired derivative of the process output 
and the actual derivative. Good control performance will be 
given by choose the following desired trajectory [3]. 



y d = a t (w -y) + cc 2 j ' (w - y)dt 



(8) 



Where a x and a 2 are design constants and w is the set point. In 
order to design a controller so that the system follows the 
trajectory defined by equation (8) as closely as possible, the 
following performance index is specified: 



} = j e 2 dt = j(y d -yy 



dt 



(9) 



i.e. minimize the error squared over a specified time horizon. In 
order to obtain y from equation (7) the chain rule must be used, 
dy dy dx 

(10) 



IV. Generic Model Control Law Design 

The choice of a controller based on Generic Model Control 
for the methane gas production control loop is based on the fact 
that the methane gas production is highly non -linear and the 
GMC allows including a non-linear model in the control 
algorithm. In this paper, we consider a biomethanation process - 
wastewater biodegradation with production of methane gas that 
takes place inside a Continuous Stirred Tank Bioreactor whose 
reduced model is presented in [1],[2]. 

These differential equations (1-5) of biomethanation 
process can be rearranged into the generic vector representation 
of the process with the following vector and scalar fields: 

State Vector: x = (X 1 S 1 ,X 2 ,S 2 ,P) 
Vector fields: /(*) = (A (x), f 2 (x), f 3 (x), / 4 (x), f s (x)) 

f 1 (x) = 0 1 -DX 1 (13) 

f 2 (x) = -Ki0! - DS 1 + DS in (14) 

Mx) = 0 2 -DX 2 (15) 

U{x) = K 2 0 1 -K 3 0 2 -DS 2 (16) 

f s (x) = K 4 0 2 - DP - Q (17) 

g vector fields: g(x) = (g x (x), g 2 (x), g 3 (x),g 4 (x), g 5 (x)) T 

The h scalar field describing the output function is simply the 
state itself, 

y = h(x) = P (18) 
These scalar and vector fields are important for the design of 
GMC law. 

With y — h(x) = P therefore, 

dh dh dh dh 

dx 1 dX t ' dx 3 dS t 

dh dh dh dh dh dh 

dx 2 dX 2 ' dx 4 dS 2 ' dx s dP 



dx dx dt 



Where 



dy 

dx 
dx 

dt 



is a row vector 



is a column vector given by eq.(l) 



dh is also a row vector and, 
dx 

dy dh ^ ^dh ^ ^ 
dt dx dx 



(11) 



Using equation (11) the performance index, the equation (10) is 
minimized when e = 0, 



U\(w — y) + a 2 j \w - y)dt - 
= 0 



dh 
dx 



dh ] 



dh dy 

— = [0,0,0,0,1] »» — = y = f 2 (x) + g 2 (x)u 

In other words, the GMC controller law is given by: 

dh 



oq (w - y) + j oc 2 ( w ~ y)dt - 
= 0 



dh 
dx 



dh ] 

/(*)+-*(*)«] 



(19) 



And solving for the manipulated variable D(t), yields, 

«1 (Pset ~ P) + j K 2 (Pset ~ P)dt - [K 4 0 2 - DP - Q] 
I 

(Pset ~ P) + j K 2 (Pset ~ P)dt - K 4 0 2 + DP + Q 



= 0 



oc 
= 0 

Then, the control law is: 



(20) 
(21) 



(12) 
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u(t) = D(t) = 



"«1 (Pset -P)-j*2 (fset ~ P ) dt 



+ ^402 - Q 



(22) 



V. Simulation Results and Discussion 



In order to examine the effectiveness of the developed 
control strategies, simulation studies were conducted using the 
model described above for the behavior of the biomethanation 
process. The model was initialized in reference to the listed 
values in [1-2] with the following (Monod-type) model for the 
reaction rates: 

S 1 -X 1 



0i = Mi 



K m - Si 



and the Haldane Kinetic model [ 1 ] 



02 = P-2 



So ' Xn 



Kmi + S 2 + S 2 /K t 



where K M1 , K M2 are Michaelis-Menten constants, n ± , \i 2 represent 
specific growth rates coefficients and K t is the inhibition 
constant. 

The performance of the generic model controller presented 
in section IV has been tested through extensive simulations by 
using the process model controller presented in section II. The 
model parameters and the nominal values of parameters of 
process considered in the simulation are: 



3.2, 
1.5, 



K t = 

n x = 0.2 hr 1 , 
K M2 = 4 g/i 



K2 
K h = 



= 16.7, K 3 = 1.035, tf 4 = 1.1935 

3, K 7 = 0.113, 

K M1 = 0.5 g/l. □ fi 2 = 0.35 h' 1 



K t = 21 g/l, 



0; S ; = 15 g/l. 



The nonlinear of biomethanation process model and control 
structure are expressed using SIMULINK block diagram, and Vi- 
rile Matlab programs. Three examples are presented to evaluate 
effectiveness of the proposed control strategies. Trial-and-error 
method is used to estimate the controller parameters oq and oc 2 . 

Figures 1 through 3 represent the controlled and 
manipulated variables response of the GMC controller to 
methane gas set point tracking and disturbance rejection. Figure 
1 illustrates the response of the GMC in tracking the methane 
gas set point, the controller performance presented fast tracking 
and no steady state error. Figure 2 illustrates the response of the 
GMC in tracking the positive change in methane gas set point at 
time 50 hr, also, the controller performance presented smooth 
behavior. Figure 3 illustrates the response of the GMC in 
tracking the negative change in methane gas set point at time 50 
hr, also, the controller performance presented smooth behavior. 
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Fig. 1 : Simulation results in tracking methane set point 
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Fig. 2: Simulation results with a change in methane set point 
from 3.0 g/l to 5.0 mg/1 
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Fig. 3: Simulation results with a change in methane set point 
from5.0g/lto3.0mg/l 

VI. Conclusion 

The GMC controller is designed and implemented to track 
the methane gas set point and to reject disturbances using a 
nonlinear model of biomethanation process. The responses of 
GMC in tracking set point are presented. It is found that the 
performance of GMC controller is more stable, more fasting, no 
overshoot. The control action has smooth behavior, but this 
behavior has sudden change at the moment of set point change 
and then return to smoothing behavior. 
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